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ABSTRACT 

The methodology for mapping ricefields in Southern 
France is developed, using 1975 LANDSAT 2 and aircraft 
data and taking into account the features of the fields. 


As part of the LANDSAT 2 programme on Agricultural Resources Investigations 
in Northern Italy and Southern France (the AGRESTE programme) (1), research has 
been carried out in the Camargue French test-site on rice crop in order to make 
an inventory of the total area of rice cultivation within the test-site. 

The study is performed at several levels : ground observations, aircraft 
and satellite imagery, after consideration of the features of French ricefielct. 

1. DESCRIPTION OF RICEFIELDS IN SOUTHERN FRANCE 

The arrangement of ricefields in Southern France depends on the ecological 
and agricultural conditions of the zone. Areas of rice cultivation consist of 
adjacent parcels, the size of which depends mainly on the local topogrfphy. The 
parcels are leveled so that water depth can be controlled to within one centi- 
meter. This implies that the parcels be narrow and that they follow contour li- 
nes. 


Given this restriction, most of the ricefields are divided into small par- 
cels (40 m, 50 m wide and 100 to 400 m long). The parcels are separated by 4 to 
10 m wide dirt roads which must permit transportation operations (sowing, chemi- 
cal treatment, aarvest..) which often take place after the fields have been 
flooded. After flooding, the fields look like a mosaic of rectangular surfaces 
with the water level at about 10 cm. 

Furthermore, several varieties are cultivated : Balilla, Euribe, Delta 
Cristal, Cigalon, Arlesienne, RB . . . and this effects the irrigation conditions 
which depend partly on the water requirements of the plants. The cover rate 
during the vegetative stage can vary also from parcel to parcel depending on 
the variety planted. 

In addi ti on ,i denti fi cati on of ricefields by their spectral properties may 
be difficult : confusion with natural water bodies at the sowing stage and 
with marshes at the vegetative stage is possible. 

2. GROUND DATA DETERMINATION 

Data were collected simultaneously with the passages of Landsat 2 for the 
vegetation season of 1975. 

The stages of development of rice are determined for the 3 most commonly 
cultivated rice varieties in the Camargue region, i.e. the Delta, Euribe, and 
Balilla 28 varieties. Table 1 summarizes the soil conditions and the growth 
stages observed during the passages of Landsat 2. The growth stages are defined 
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according to international standard (2). 


At the same time climatic measurements were recorded in an experimental 
ricefield giving air temperature, wind (direction, speed) soil temperature, 
sun, relative air moisture and rainfall recordings (3). 

In conjunction with ground observations and measurements, Landsat 2 data 
processing plan can be outlined (table 2). 

3. STUDY OF THE AERIAL IMAGERY 

During a B 17 flight over the test site at 7000 m and 1500 m, June 20, 197/, 
IRC and IR black and white images, as well as Daedalus Scanner recordings were 
obtained. 

Photo interp re tation of 7000 m IRC images provides the location of several 
samples of ricefields and other field categories in the area : vineyards, wheat- 
fields, marshes and urban areas. On June 20, the ricefields are under water and 
different varieties exhibit different growth stages : from emergence to begin- 
ning of tillering. On IRC images ricefields appear as blocks of mixed blue and 
pink tone. 

Daedalus images, especially those from the 2 channels closest to the MSS 5 
and MSS 7 bands of LANDSAl : channel 7 : 650 - 700 nm, channel 10 : 900 
1100 nm are studied. On channel 7, the interpretation presents some ambiguity 
between ricefields and wheatfields while in channel 10, ricefields are confused 
with temporarily flooded plots (vineyards, market gardens). 

In conjunction with photointerpretation and knowledge of the topography, the 
pedological and hydorlogical considerations of the zone, the test site which 
includes the totality of ri cecul ti vati on in France is devided into 4 rectangular 
areas or strates, suitable for computer aided classification. 

Figure 1 shows the test site (43° 24' to 43° 45' North, 4° 25' to 4° 50' 
East) divided into 4 subsets : 


1. 

Arles Area : 

20 

km 

X 

20 

km 

2. 

Aigues Mortes : 

20 

km 

X 

11 

km 

3. 

Ste-Marie de la 

Mer: 11 

km 

X 

11 

km 

4. 

Port St-Louis : 

15 

km 

X 

20 

km 


4. METHODOLOGY FOR DATA PROCESSING 

On the 8 LANDSAT images expected, 4 were available ; July 6 (2165 - 09511), 
July 23 (2182 - 09451), August 11 (2201 - 09502), October 3 (2254 - 09440). 
Landsat and aircraft data are processed first for the study of ricefield res- 
ponse. Then Landsat data are processed for ricefield classification. 

4.1. STUDY OF RICEFIELD RESPONSES 

Photo-i nterpretatoon of aerial images serves as a guide to the study of the 
homogeneity of ricefields responses. 

An area is defined as homogeneous when the gradient at all points is zero 
or less than a threshold determined by the measurement noise. In addition to 
this concept of instrumental homogeneity, there is also the concept of the 
homogeneity of the ground object. The latter is related to the resolution 
problem. With Landsat resolution, such agricultural terrain as vineyards or 
orchards with a structure of alternating vine and soil or trees and soil appear 
as stable composites and the samples taken within these types of terrain have 
homogeneous responses. 

4.1.1. FROM AERIAL IMAGES : to study the homogeneity of ricefields, 1500 m 
IR black and white images are digitized by means of the Joyce-Loebl Scandig 
microdensitometer, delivering an 8 bit coded optical density for every pixel 
of about 2 m X 2 m. 
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Figure 2 shows a part of an original image including ricefields, wheafields 
and meadows. Several samples of each field category are located and studied by 
means of their histogram. As expected, the histogram of ricefield samples 
exhibits larger dispersion than other fields. Table 3 shows some statistics 
computed on ricefield, wheatfield and meadow samples, each of which corresponds 
to an area of 40 pixels x 30 pixels (80 m x 60 m). The case of Landsat pixel can 
be then simulated by taking the mean value of each sample. The non uniformity of 
ricefield responses is then demonstrated by the wide range of mean values of 
ricefield samples in table 3. 

4.1.2. FROM LANDSAT IMAGES : Several samples are extracted from Landsat 2 
images of July 6, July 23 and August 11 : cereals, vineyards, meadows, rice- 
fields. Figure 3 shows the Landsat 2 image enlargement of the test site (on 
July 23, 75) and table 4 compares statistics computed from 3 samples of rice- 
fields and 3 samples of wheatfields of the same size (about 60 pixels) on J ly 
23 and August 11. Figure 4 shows examples of ricefield and wheatfield histograms 

The wide range of ricefield responses in June and July can be explained by 
the different growth stages of the rice plants in different parcels : beginning 
of the tillering, beginning of the stem elongation, boots just visible. The 
water layer appears through the vegetation and the response of one parcel can 
differ from that of its neighbour. In August, ricefields are at the beginning of 
the heading, the fields are then well cowered by the vegetation and the respon- 
ses are more uniform. 

In figure 5, several field categories are studied by diagramming MSS 5, 

MSS 7 sample responses from the July 6, July 23 and August 11 images. Responses 
of ricefield appear in July images as a cluster which can be separated from the 
others, meanwhile for August date, ricefields and vineyards are merged in one 
cl uster. 


In the vegetative stage of the rice plant, there is a correlation between 
responses of vegetation and water, and the correlation coefficient 
oXi y. 

r = is greater than that of other field categories, where a„ and o„ are 

Ox-ay y X y 

standard deviations of MSS 5, MSS 7 responses. As an example, correlations 


coefficient computed from July 6 

are the 

following : 



°MSS 5 

°MSS 7 

r 

ri ce 

4.37 

6.44 

- 0.117 

vine 

4.48 

3.57 

- 0.012 

wheat 

7.57 

5. 15 

- 0.005 


For the July images, the fastest method for distinguishing ricefields can 
consist of choosing a linear classifier in the MSS 5 - MSS 7 plane. 

4.2 DATA REDUCTION 

To reduce the dimension of mul ti temporal Landsat data, Sebastien discrimina- 
tion criterium (4) is applied on samples of the main field categories : rice- 
field, wheatfield, vineyard extracted from July 6, July 2? and August 11 
Landsat images. 


The samples are divided arbitrarily into 


(x. Y,^) 


1/n 


( (X - E(Y|^) )^ S'^ (X - 


with Y, 


set of samples (k = 1, 2, 3) 
covariance matrix of Y|^ 
number of channels 


2 groups. The 
E(Y|^) ) + n ) 


( 


discriminant function 


E : expectation function 

computed on the first group, is then applied to the second for calculation of 
the "wel 1 - cl assed" sample percentage. Combinations of 1, 2, 3, 4, 5... channels 
among the 12 are computed in order to determine the best one for discriminating 

1485 


ricefield, wheatfield and vineyard. 


Table 5 results calculation of "we 1 I -cl assed" percentage with data 
from different date and band and shows in the 3 cases the best combinations 
of 2, 4, 5 channels among the 12. 

The August image appears the best for ricefield recongnition : the rice- 
field response is then well separated from that of wheatfield (harvested) and 
vineyard. 90 % of ricefield can be computed from only MSS 5, MSS 7 images. 
Unfortunately, for inventory purpose, the August image cannot be used without 
combination with an other date data, the test site being partly cloud covered. 

4.3. CLASSIFICATION METHODS 

Unsupervised and supervised methods have been applied to a limited test- 
site for evaluating the efficiency of the methods compared with ground determin 
nations. 

4.3.1. UNSUPERVISEO METHOD : the unsupervised method evolved is a mobile 
center clustering technique called “Nuees Oynamiques" by E. Oiday (5). The 
principle of the method consists of settling clusters around centers Gj chosen 
randomly on the image. In the next step the process is iterated, i.e. 
new centers , g 3 ... G^ must be built and stopped when the centers ramain 

fixed, i.e. when G" = G’J'*’! VC. 

1 ? 

THe method is justified by the fact that the passage from G^ to G^., Gc. . . 
reduce the intra-class variance of E^-. 

To reduce the number of iterations, the area under study is divided into 
subareas, each of which_is a square of 60x 60 pixels. 


The process is applied first to Z... The iteration number is 
limited to 5. If after 5 iterations, the convergence of the 
process is not reached, an other random choise is made. When 
the convergence on Z., is obtained, zone Z,- is computed with 
the last set of centers of Zjj. 


'll 

'12 

^13 

'21 

'22 

'23 

'31 

^32 

'33 


In the case of Landsat imagery, every sub-area corresponds to a zone of 
about 4.8 km x 3.6 km. Geographically, in such an agricultural zone, 2 conti- 
gous subareas are in most case included in the same ecological region. This 
means that the same types of vegetation or land use will very probably be 
found in 2 contiguous subareas. So, in practice, one iteration is necessary for 
computing 2^2 with the last centers of Z^j. If however the "landscape" changes 

suddenly, then other iterations are needed, or a choice of random initial center 
must be made. The contiguity concept can be applied in both directions, i.e. 
lines and col urn ns. Thus for Zgj, the initial centers came from Z ^ . 

4.3.2. SUPERVISED METHODS : several methods are applied, but parametric 
method such the maximum likelihood ratio is not utilized by the fact that the 
normal hypothesis for ricefield responses is not proved and the a priori occu- 
rence probability of classes are unknown. 

- Barycentric method : using weighed euclidean distance as discriminant 

function : ^ 

d (x^, X 2 ) = k-j ■ ^^2i^ 

: weighting factor for channel i. The principle consists to chosseK.,- 

to get a maximum of well classed prototypes. The process starts with 
112 

= 2 • where is the total variance of prototype in channel i. 

°i i 

- Quadratic method : using Sebastien (4) discriminant function. 

- Fix-Hodges method (6) : using Mahalanobis distance and generalizing the 
K-Neares t-Nei ghbourg rule. - 

- Chi -Square-method : using XT distance as classifier : 
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2 


d ( Xj , 


■± 

i =1 



^2i 

n 

2 X2i 


4.3.3. COMPARATIVE STUDY : to illustrate the efficiency of the methods, dif- 
ferent algorithms are performed on August 11 LANDSAT data, pointed out as the 
best date for ricefield inventory, according to results of the data reduction 
study. 

Figure 6 shows a 6.7 x 9.5 km area of the test site (120 pixels x 120 pixels) 
with the location of some field categories : ricefields in dark grey, vineyards 
in medium grey and wheatfields in light grey. Different results are studied : 

- Barycentric and quadratic methods : the reference ricefields are mapped 

{ f i gure 7) wi th 3 of difference between computed surface and owner estimation. 

Vineyards are subdi vised into 2 classes, relating to their soil conditions. 

- Fix-Hodges method : ricefields are mapped with 5 % of divergence. Vine- 
yard is classified as one class, including meadows, market garden and bare soil. 
For ricefield, the method can fit to their non normal spectral responses. 
However, it requires a very good sample selection which is difficult to per- 
form on such parcelled agricultural site. 

- Unsupervided methods : results are obtained with clustering techniques 
using euclidean, Mahalanobis and distance. The ricefields and vineyards 
are partly merged as expected, according to diagram 5C, while wheatfields are 
divided into 2 classes which do not correspond to difference in cultural 
condi ti ons . 

Figure 8 shows an example of unsupervised method applied on the site, using 
Mahalanobis distance. 

Table 6 compares areas in hectares of some most important fields computed 
by different techniques with owner estimations. 

The computer CP time for supervised methods is 0.5 sec for 60 x 60 pixels 
and 5 classes. However, for Fix-Hodges discrimination, the training phase for 
750 samples is about 5 sec. For clustering techniques, computer CP time is 
function of the iteration number, each iteration taking about 1 sec. In prac- 
tice, 3 iterations are required for the first subarea and one iteration for 
the followings then for 120 x ,120 pixels, about 10 sec are required. 

4.4. TEXTURE ANALYSIS : 

Texture analysis is attempted on 1500 m IR B&W imagery by edge detection and 
study of the texture within detected fields. 

4.4.1. EDGE DETECTION : the algorithm uses discrete operators to detect 
edges of several fields and marshes in the image. The gradient of grey level 
is replaced following Eber 1 ei n-Bes zk a ' s method by (7) 

Gj^(i,j) - DIF (i,j) , where ^ x 

Gx(i»j) = ( (a-x (b-x)2 ) and d c 

DIF (i,j) =x-x, x=(a+b+c+d+ x)/5 

The edge detection is computed by an algorithm starting on the maximal 
gradient point on the image and following the line of minimal slope. The chain 
is stopped when either a low threshold, or a point belonging to a previous 
computed chain is reached. 

Figure 9 shows result of edge detection technique computed on part of the 
aerial image (figure 2). 

4.4.2. TEXTURE PROPERTY : by analyzing the frequency distribution of the 
values of the field local porperty, one can obtain invariant properties of the 
field. Statistics of local porperties - called textural properties - can be 
served as a decision rule. In table 3 and 4, parcels of early growth stage 

1487 



riceplants (June, July) are characateri zed by t a wide range of their responses 
with regard to the other parcels. This property is used to recognize ricefields. 

4.5. EVALUATION OF UNPRODUCTIVE SURFACE : 

With the resolution of Landsat, ricefields are mapped including rc s and 
dikes (4 - 5 m width). For a statistical evaluation of the "unproductive" 
surface, several ricefields are studied from 1500 m aircraft image (June 20) by 
means of histograms. In most case, the road responses differ from that of 
flooded rice ields and a density slicing provides ratio of road area on the 
total. The percentage averaged over 9 representative ricefields is 17,2 %. 


4.6. CONCLUSION 

For the best use of Landsat and aerial data for ricefield mapping purpose, 
the conclusion drawn from the 1975 data on the French test-site are as follows : 


June 20 : - recongnition of ricefield by edge detection and texture 

analysis (aerial image) 

- evalua ion of unproductive surface within a block of ricefield 
(aerial i mage) 

July 6 : - cl assi fication by linear di s cr imi na t i on (LANDSAT) 

July 23 : (confusion with marshes is possible) 

August 11 : - classification by barycentric or quadratic distance 


Jul . 6- Aug . 1 1 : 


Jul.6-Jul.23-, 
August 11 


classification using barycentric distance on MSS 7 - MSS 5 
bands of the 2 dates 

classification using barycentric distance on MSS 7 - MSS 5 
bands of August 11, MSS 7 band of July 6, MSS 5 band of 
the July 23 images. 


The inventory method using barycentric distance on the last data combination 
gives about 12 ODD hectares within the test-site, while the value given by the 
agricultural statistics is 10 500 hectares. Figure 10 shows result obtained on 
the s ub- tes t- s i te N° 1 (Arles area), the darkest sign corresponding to ricefield. 
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date of 
passage 

Delta 

growth stage 
E uri be 

Balilla 28 

soil condi ti on 

April 7 

- 

- 

- 

wet 


Apri 1 25 

- 

- 

- 

dry 


May 13 

Emergence 

Erne rgence 

Erne rgence 

under 

water 

May 31 

Emergence 

Emergence 

Emergence 

M 

II 

June 18 

Beginning of 
ti 1 leri ng 

Emergence 

Emergence 

II 

II 

July 6 

Beg. of stem 
el ongati on 

Beg. of 
tillering 

Beg. of 
ti 1 leri ng 

II 

II 

July 24 

Boots just 
visible 

Beg. of stem 
e 1 ongati on 

Beg. of stem 
e 1 onga ti on 

II 

** 

Aug. 11 

Mi d 

f 1 oweri ng 

Beg. of 
heading 

Boots just 
visible 

II 

II 

Aug. 29 

Early milk 
stage 

Mi d 

f 1 owe ri ng 

Mid 

f 1 oweri ng 

II 

II 

Sept. 16 

Early dough 
s tage 

Early milk 
stage 

Early milk 
s tage 

II 

II 

Oct. 4 

Caryopsi s 
loosening at 
day ti me 

Caryopsi s 
hard 

Early dough 
stage 

Partly 
1 evel 

at zero 

Oct. 22 


C a ry 0 p s i s 
loosening at 
day ti me 

Caryopsis hard 



Nov. 9 


Ha rves ted 

Harvested 

- 



TABLE 1. GROWTH STAGES OF RICEPLANTS AND 
SOIL CONDITIONS during the pas- 
sages of LANDSAT 2, 1975. 



es 

Object! ve 

remarks 

May 

13 

Ricefield inventory 

Result includes natural water bodies 

May 

31 

II 

II 

Perturbation : rain gauge reading = 11 mm 

June 18 

II 

II 

Several r.icefields were dried out for wee<l 
killer treatment 

July 6 

II 

II 

Result includes ricefields ans marshes 

July 24 
August 11 

Sept. 16 
October 4 

M II 

Ricefield inventory 
and disease dete-o- 
tion 

II II 

Mapping of rice 
varieties 

Result includes ricefields, marshes and 
temporarily flooded plots (vineyards, 
wheatfields after harvest...) 


TABLE 2. 
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SAMPLE 

MEAN VALUE 

STANDARD DEVIATION 

Ri cefi e 1 ds 

19.32 

9.39 


39.98 

10.23 


18. 83 

5.57 


32.51 

11.45 


15.95 

5.29 

Wheatf i el ds 

20.07 

2.64 


17.89 

3.87 


20.48 

4.60 

Meadows 

28.48 

3.29 


25.52 

4.41 


23.43 

4.06 


TABLE 3. STATISTICS OF SOME FIELD RESPONSES 
June 30. 1975. 


Date 

samp 1 es 

MSS 

mean value 

5 

standard 

deviation 

MS 

mean value 

5 7 

sta ndard 
deviation 

July 23 

Ri ce 

28.86 

5.46 

31.79 

3.95 

1975 


23.58 

6.31 

29.08 

5.09 



33.50 

4.68 

31.08 

4. 12 


wheat 

43.04 

2.10 

23.03 

1.16 



46.71 

2. 16 

27.75 

1.03 



41.71 

1.95 

24.20 

1.38 

Aug. 08 

Ri ce 

22.31 

1.27 

36.18 

3.26 

1975 


19.95 

1.30 

35.75 

2.45 



18.29 

1.73 

37.48 

2.54 


wheat 

38.92 

4.01 

46.86 

3.92 



40.50 

3. 13 

51.21 

3.43 



31.21 

2.92 

37.85 

5.80 


TABLE 4. STATISTICS OF SOME FIELDS LANDSAT 
RESPONSES 
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MSS 7 

, MSS 

Aug . 1 1 

MSS 

MSS 

MSS 

1 , MSS 5 Aug. 11 
July 6 
■ August 11 

MSS 

MSS 

MSS 

MSS 

7, MSS 5 Aug. 11 
7 July 6 
August 11 
July 23 

CQjii^puted 
ori gi^a 

ri ce 

wheat 

vineyard 

ri ce 

wheat 

vineyard 

ri ce 

wheat 

vineyard 

ri ce 

209 

11 

12 

214 

5 

13 

214 

5 

13 

wheat 

3 

172 

33 

12 

176 

20 

7 

183 

18 

vineyard 

12 

7. 

79 

3 

9 

86 

3 

96 

89 


well classed percentage : 
85,50 * 

well classed percentage 
for rice : 90 % 

well classed percentage : 
88.48 % 

well classed percentage 
for rice : 92 % 

well classed percentage ; 

90.33 % 

well classed percentage 
for rice : 92 * 


TABLE 5 



Les Pebrieres 

Les Huits Clos 

Mas du Tort 

Mas de Vert 

Dorn, de la Reine 


Rice 

Ri ce 

Ri ce 

Vineyard 

wheat 

Owner estimation 

43 

189 

95 

84 

111 

Barycentri c 

45 

177 

82 

82 

113 

Quadratic 

47 

195 

86 

73 

115 

Fix-Hodge 

46 

198 

71 


117 

Clustering- 
euclidean d. 

46 

190 

94 

72 

52 + 67 

Clustering- 
Maha lanobis 

42 

157 

66 

73 

88 

Clustering- 
X2 d. 

45 

157 

87 

77 

53 + 67 

Cl ustering- 
4 classes 

44 

198 

87 

63 

102 

Clustering- 
6 classes 

40 

106 ^ 
113 

17^ 

58 

41 + 62 


TABLE 6. COMPARISON OF FIELD AREA COMPUTED FROM 
LANDSAT DATA WITH OWNER ESTIMATION. 
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FIGURE 1 . CAMARGUE FRENCH TEST SITE 



BS! 

FIGURE 3. LANDSAT 

IMAGE ENLARGEMENT 


MSS 7 - 

July 23 - 

FIGURE 2. PART OF 1500 
IMAGE. 

( 1 ) : block of 

m AERIAL 
ricefield. 

Scale ^ 

1/750 000 . 





MSS5 


a : July 6 


b : July 23 c : August 11 

t 

FIGURE 5. MSS 5 - MSS 7 response diagrams/ 

rice, ... wheat, . _ ._vi neyard. 

marshes 
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FIGURE 6. GEOGRAPHICAL MAP. 

Part of the test-site. 


FIGURE 8. CLASSIFICATION RESULT : 
unsupervised method. 


FIGURE 7. CLASSIFICATION RESULT : 
supervised method, 

barycentric distance, dark grey;rice. 



FIGURE 9. RESULT OF EDGE DETECTION 
TECHNIQUE applied on 
figure 2. 
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FIGURE 10. CLASSIFICATION RESULT : barycentric distance. 

July and August data. Dark gray: ricefield. 
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